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Transport properties of arrays of metallic quantum dots are governed by the distance-dependent 
exchange coupling between the dots. It is shown that the effective value of the exchange coupling, 
as measured by the charging energy per dot, depends monotonically on the size of the array. The 
effect saturates for hexagonal arrays of over 7 5 unit cells. The discussion uses a multi-stage block 
renormalization group approach applied to the Hubbard Hamiltonian. A first order phase transition 
occurs upon compression of the lattice and the size dependence is qualitatively different for the two 
phases. 
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Two-dimensional arrays with tailored electronic prop- 
erties are generating much current experimental and the- 
oretical interest [1-6]. This was made possible by the de- 
velopment of synthetic methods for the preparation of, 
so-called, quantum dots (QDs). There is much current 
interest in the ability to engineer nanoscale electronic de- 
vices, including quantum computers from quantum dots. 
Quantum dots of nearly identical sizes self assemble into a 
planar array. (The dots are passivate against collapse by 
coating them with organic ligands). For, e.g., Ag nan- 
odots, the packing is hexagonal. The lower lying elec- 
tronic states of an isolated dot are discrete being de- 
termined by the confining potential (and therefore the 
size) of the dot. Because of their larger size (100-1000 
atoms each) it takes only a relatively low energy to add 
another electron to a dot, as revealed by scanning tun- 
neling microscopy [7]. This energy is much lower than 
the corresponding energy for ordinary atoms and most 
molecules. It follows that when dots are close enough to 
be exchange coupled, which is the case in an array, the 
charging energy can be quite low. We here propose and 
implement a computational method that allows the con- 
tributions of such ionic configurations even for extended 
arrays. The technical problem is that the Coulombic re- 
pulsion between two electrons (of opposite spins) that 
occupy the same dot cannot be described in a one elec- 
tron approximation. It requires allowing for correlation 
of electrons. Most methods that explicitly include corre- 
lation effects scale as some high power of the number of 
atoms (here, dots) and are computationally intractable. 
For example, a hexagonal array of only 19 dots, 3 dots 
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FIG. 1: Schematic diagram of the triangular lattice with 
hexagonal blocks. Only two neighboring blocks p and p' are 
drawn here. The dotted lines represent the interblock inter- 
actions and solid line intrablock ones. 



per side, has already 2,891,056,160 low electronic con- 
figurations. So earlier [8] exact computations including 
charging energy were limited to a hexagonal array of only 
7 dots, 2 dots per side. Yet current measurements of both 
static [9] and transport [10] properties use arrays of at 
least 100 dots per side. The simplest Hamiltonian that 
includes both the Coulombic (or charging energy) effects 
and the exchange coupling is the Hubbard model [11]. 
This model can be solved exactly for a one-dimensional 
chain, but for a two-dimensional array it is, so far, ana- 
lytically intractable. In the absence of a closed solution, 
various methods have been developed [12]. Renormal- 
ization group (RG) methods are receiving increasing at- 
tention because of their non-perturbative nature, that 
allows application to the intermediate-to-strong coupling 
regime. This Letter reports the applications of a real- 
space block renormalization group (BRG) method [13- 
15] on a 2-dimensional triangular lattice with hexagonal 
blocks as shown in Fig. 1. This enables us to explore the 
evolution of size effects. Specifically, we demonstrate how 
the effective exchange coupling between adjacent dots is 
modulated by the size of the array. In other words, we ar- 
gue that for finite arrays of quantum dots, their response 
will be size dependent and we determine the scaling. 
The essential physics of the collective low energy elec- 
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tronic structure of an array of metallic quantum dots 
has the following ingredients. Each dot is characterized 
by the energy \i needed to remove its highermost (=va- 
lence) electron. Adjacent dots are exchange coupled with 
a strength, t, that decreases exponentially at larger dis- 
tances but saturates at very close packing. Transport 
processes are facilitated by the coupling and are therefore 
measured in the intermediate to strong coupling range, 
U jt. When t is large it can overwhelm the role of the 
charging energy U and it is then possible to treat the 
lattice by a tight binding approximation that can be eas- 
ily applied to larger size arrays. However, over much 
of the experimentally available compression range it is 
the case that t < U. Our letter specifically addresses 
this intermediate coupling range where the correlation of 
electrons needs to be explicitly accounted for. 
The Hubbard Hamiltonian is written as[16] 

H = -t ^2 4a c 3« + U Y1 n *t"a _ ^X! + n a) 

<i,j>,a i i 

(i) 

c| CT (c i(T )creates(annihilates) an electron with spin tr in the 
valence (=Wannier) orbital of the dot located at site i; 
the corresponding number operator is rii a = c\ a Ci a . The 
angular bracket < ... > on the first sum in Eq. (1) in- 
dicates that summation is restricted to nearest-neighbor 
dots. This model Hamiltonian allows only one orbital per 
dot. That orbital can be empty or it can accommodate 
one or two electrons. 

The essence of the BRG method is to map the above 
many-particle Hamiltonian on a lattice to a new one with 
fewer degrees of freedom and with the same low-lying en- 
ergy levels [17]. Then the mapping is repeated leading 
to a final Hamiltonian of a seven site hexagonal array for 
which we obtain an exact numerical solution. The pro- 
cedure can be summarized into three steps: First divide 
the TV-site lattice into appropriate n s -site blocks labeled 
by p(p=l,2,..., N/n s ) and separate the Hamiltonian H 
into intrablock part Hb and interblock Hjb ■ 

H = H B +H IB =Y,H P + Y, V p-p' ( 2 ) 

p {p,p') 

where H p is the Hamiltonian (1) for a given block and 
the interblock p,p' coupling is defined in (3) below. 

The second step is to solve H p exactly for the eigen- 
values E Pi and eigenfunctions $> p i(i = 1,2, ...,4™ s ). Then 
the eigenfunctions of Hb are constructed by direct mul- 
tiplication of The last step is to treat each block 
as one site on a new lattice and the correlations between 
blocks as hopping interactions. 

The original Hilbert space has four states per site. If we 
are only concerned with lower lying states of the system 
as when studying the metal-insulator-transition [18], it is 
not necessary to keep all the states for a block. 

To make the new Hamiltonian tractable the reduction 
in size should not be accompanied by a proliferation of 



new couplings. Then one can use an iteration procedure 
to solve the model. To achieve this it is necessary to keep 
only 4 states in step 2. Their energies are £^=1,2,3,4). 
In order to avoid proliferation of additional couplings in 
the new Hamiltonian the four states kept from the block 
cannot be arbitrarily chosen. Some definite conditions as 
discussed in Ref. [16] must be satisfied. For example, the 
states must belong to the same irreducible representation 
of Cqv symmetry group of the lattice. In particular, in or- 
der to copy the intrasite structure of the old Hamiltonian, 
one must have that E% = E±. Furthermore, particle-hole 
symmetry of a half- filled lattice requires that E\ = E 2 - 
Further restrictions follow from the need to make extra 
couplings vanish. Operators in the truncated basis are 
denoted by a prime so that the interblock coupling of 
Eq. (2) is 

V pp , =v\HY J 4^v'o, (3) 

(7 

where v represents the number of couplings between 
neighboring blocks. The coupling strength for the bor- 
der sites of a block by A and the renormalization group 
equation for the coupling strength is 

t' = v\ 2 t, (4) 

The other renormalization relations are 

U' = 2{E 1 -E 2 ), K' = (E 1 +E 2 )/2. (5) 

K sets the zero of energy with the initial value of K = 
-17/4. For the half-filled lattice p = U/2. Below we 
will compute the energy gap for adding or removing an 
electron from a half-full lattice. The gap is independent 
of the choice for p. The results of the procedure are 
applied to a finite array of 7™ sites in Fig. 2. The coupling 
constants for the initial array are U = Uo and t = 1 so 
that the abscissa spans the intermediate coupling range 
Uq > to. The ordinate gives the renormalized coupling 
constants, U'/t', on a logarithmic scale, for various values 
of the size, n, of the original lattice (n = 1 being shown 
for reference because it is the line U'/t' = Uo/to). There 
are two features to note. One is that there is a phase 
transition at Uo/to = 12.5. This is quite near the recent 
result Uo/to = 12.07 obtained with exact diagonalization 
method [19] . The other feature of the results is the one we 
want to focus on. It is the variation of the renormalized 
coupling constants with size of the array. The reason for 
drawing attention to this scaling is as follows. 

The low lying electronic states of the large array are 
computed by an exact diagonalization [8] of the Hubbard- 
like Hamiltonian of a seven site array, the smallest hexag- 
onal structure, with the renormalized coupling constants. 
In the (albeit, approximate) real space group renormal- 
ization procedure that we use the only way that the size 
of the original array comes in is in the value of the renor- 
malized coupling constants. Arrays of all sizes can there- 
fore be made to fall on a common universal plot where 
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FIG. 2: The renormalized coupling strength, U' /t', logarith- 
mic scale, starting from a lattice of n hexagons of seven dots 
each, vs. the initial value of the charging energy Uo/to, (mea- 
sured in units of the exchange coupling to so that the renor- 
malization equations begin with to = 1)- The solid line is 
that for a single hexagon for which U'/t' — Uo/to- The re- 
sults show that there is a phase transition at Uo/to = 12.5. 
Note how the size effect is different at the two sides of the 
transition and how the role of the exchange coupling is in- 
creasingly important for compressed (Uo/to smaller) lattices. 



they differ only in the input values of the coupling con- 
stants. To be sure, the actual value of the coupling 
constant U'/t' depends on n, and quite systematically 
and dramatically so, Fig. 3, but for given initial cou- 
pling constant, Uo/to, there is no other dependence. In 
our earlier work we have argued, on the basis of pertur- 
bation theory and also exact computations that in the 
regime t > U the effects of the charging energy can be 
neglected, particularly so in application to larger arrays 
[6, 9, 20]. Fig. 3 indicates that this should be the case 
anywhere to the left of the phase transition, where U'/t' 
decreases strongly with increasing n. As discussed be- 
low, the low lying electronic states are determined by the 
value of U'/t' and therefore, in the intermediate to strong 
coupling regime the importance of the dot-dot exchange 
coupling increases with the size of the lattice. Conversely 
so in the weak coupling regime. To emphasize the uni- 
versal scaling we discuss the gap A between the highest 
occupied and lowest empty state of the entire array. The 
closing of this gap is the signature of the Mott insulator 
to metal transition. We define the gap as usual by the 
energy cost difference between adding or removing one 
electron from the half-filled array. This gap is the limit- 
ing value of U' after a large number of iterations but it 
can also be computed for any finite value of lattice size 
n. As is to be expected, such a plot of the scaled gap, 
A n /t' vs. Uq looks remarkably similar to Fig. 2 that 
showed the scaled coupling constant, U'/t' vs. Uo- In 
particular, such a plot establishes the phase transition 
at Uo/to = 12.5. Here we present it as a universal plot, 
Fig. 4, showing the scaled gap vs. the scaled coupling 
constant. The results for the gap for different values of 
n (identified in the legend) all fall on the same curve. 



FIG. 3: The renormalized coupling strength, U'/t', logarith- 
mic scale, for a lattice of n hexagons of seven dots each, vs. 
n for different values of the charging energy Uo of the full lat- 
tice. Note the qualitatively different trends in the weak and 
the intermediate coupling regimes. 




FIG. 4: The gap (^difference between the energy needed to 
add or to remove an electron from a half-filled lattice) for a 
lattice of n hexagons, logarithmic scale, vs. the renormalized 
coupling strength. The discrete points are computed for dif- 
ferent values of n and charging energy Uo of the full lattice. 
The line is the exactly computed value of the gap for a single 
hexagon, for which U'/t' — Uo/to- The kink is the remnant 
of the phase transition for n — > oo. 

The line going through the individual points is not a fit 
but is the results of exact computations of the gap for 
a hexagonal array of seven dots, n — 1. Larger arrays 
are renormalized first and then the results are put in as 
discrete points, corresponding to the series of values of 
Uo shown in Figs. 2 and 3. The kink in the plot is the 
memory of the phase transition for n — > oo. The points 
for n > 1 are all for the same common set of values of 
Uq and their spread on this plot is another manifestation 
of the size effect. Specifically note how points for larger 
values of n congregate toward the two extremes of the 
plot. 

Fig. 4. shows that the value of Uo/to f° r which the 
transition occurs depends on the size of the lattice. The 
value (Uo/to — 1.5) is quite low for n = 1 and increases 
with n saturating at 12.5 for n — » oo. The physics of this 
increase is that originally made by Mott [21] in identify- 
ing the insulator to metal transition. It is a competition 
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between localization, measured by U , and derealization, 
measured by t. Other things being equal, the bigger the 
lattice the more an electron can delocalize. From the 
point of view of renormalization, at each iteration we 
lump seven dots into one. Such a bigger dot has a lower 
charging energy and furthermore a renormalized dot is 
more strongly coupled to its neighbors, cf. Eq.(4). Both 
effects reduce the critical value of U~o/t . It should be 
mentioned here that in order to increase the accuracy, the 
hexagonal blocks instead of the most-used triangular ones 
are utilized. In fact, for triangular blocks, similar scalings 
are also found. There is no qualitative difference, which 
gives us much confidence upon the consistence of the re- 
sults between the RG procedures with difference kind of 
blocks. Nanoscale devices are finite arrays of quantum 
dots. The results of the renormalization group approach 
suggest that such arrays can be usefully treated and dis- 
cussed by renormalizing the effective coupling strength 
down to the unit cell of the array. The scaling with size 
is predicted to vary depending on the coupling regime. 
The dot-dot exchange coupling becomes of lesser impor- 
tance for larger arrays that are weakly coupled (or have 
a higher charging energy) and it is of increasing effective 
strength for intermediate to strong coupling. The avail- 
able computational and experimental evidence [6, 9, 22] 
is that the coupling regimes and the resulting phase dia- 
gram are far richer than the single phase transition that 
we here discussed. This is for two reasons. One is that 
there are additional control variables (e.g., external elec- 
trical field [9], temperature [6]) that we have not intro- 
duced here. The other is the inherent size fluctuations 
of quantum dots [22] that can induce a transition to a 
domain localized phase [23]. It corresponds to charge 
transfer between not near neighboring dots and is sim- 
ilar to the super-exchange coupling in molecular charge 
transfer. This transition where the electronic states are 
not localized but are neither extensively spread out seems 
to have been experimentally detected [9] . We plan to ex- 
tend the RGB procedure so as to allow us to treat the 
role of disorder. Preliminary results already verify that 
this second transition can be seen by renormalization of 
larger lattices. Equally we seek additional experimen- 
tal signatures of the variation of the effective coupling 
strength with the size of the array. In conclusion we note 
that a single quantum dot is made up of tens to hundreds 
of atoms (or molecules). Yet its very low lying electronic 
states can be simply and usefully described by the idea 
of a quantum confinement. Here we tried to do the same 
for a finite array of quantum dots. A larger number of 
hexagons was scaled down to a single hexagon with an 
effective coupling. This scaling is achieved by collapsing 
seven dots into one at each stage of the iteration. 
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